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Abstract: The Morato's digger toad {Proceratophrys moratoi) inhabits Brazilian moist 
savannas and is critically endangered due to its very limited geographic distribution, 
reduced number of isolated populations, and evidence of population decline and local 
extinctions. With the objective of providing tools for the genetic study of the species, 
22 polymorphic microsatellite loci were isolated and screened using DNA extracted from 
samples of oral mucosa cells obtained from 113 individuals representing five remnant 
P. moratoi populations in the Brazilian state of Sao Paulo. These markers presented 
2-18 alleles per locus, polymorphism information content (PIC) of 0.02-0.87, observed 
heterozygosity of 0.02-0.96 and expected heterozygosity of 0.02-0.87. Three of the loci 
deviated significantly from Hardy-Weinberg equilibrium in one of the populations, 
possibly due to the presence of null alleles. Significant linkage disequilibrium was also 
detected between three pairs of loci. The molecular markers developed in this study were 
able to discriminate each of the individuals sampled (identity analysis). This means that 
they will be extremely useful for future genetic studies applied to the conservation of 
P. moratoi, providing a baseline for estimating the levels of genetic diversity, pedigrees, 
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inbreeding, and population structure, which will be essential for the development of 
effective genetic management programs. 

Keywords: Proceratophrys moratoi; endangered species; microsatellite; population 
genetics; Morato's digger toad; conservation genetics 



1. Introduction 

The class Amphibia has experienced a major global decline in recent decades, becoming more 
endangered than birds and mammals, due to a combination of factors [1]. Habitat destruction, climate 
change and infectious diseases are considered to be the primary cause of the decline of this group [2,3]. 
Ongoing anthropogenic impacts have contributed to the increasing deterioration of landscapes, which 
not only modifies aquatic and terrestrial habitats, but also reduces their connectivity, which are all 
factors that may affect amphibian populations adversely [4,5]. 

Proceratophrys moratoi is a digger toad of small size, typically with a snout-vent length of no more 
than 35 mm, which is endemic to the Cerrado savanna of the Brazilian state of Sao Paulo [6]. The 
species is found in campo sujo habitats (grassland dotted with small shrubs), invariably near the 
gallery forests associated with the headwaters of streams [6,7]. In Sao Paulo, the Cerrado biome has 
been modified intensively in recent decades, primarily for the planting of commercial crops such as 
sugarcane, but also for cattle ranching and urban development [8]. Currently, only about 6% of the 
original cover remains [9], which has drastically reduced the availability of potential habitat for the 
endemic P. moratoi. Due to its very restricted geographic distribution and the evidence of population 
decline and local extinctions [7,10], the species is currently listed as critically endangered by the 
International Union for Conservation of Nature [11], and is included in the official lists of endangered 
species of Brazil [12] and Sao Paulo [13]. 

The available genetic studies of P. moratoi include molecular analyses of mitochondrial and nuclear 
genes [14] and cytogenetics [15]. However, no population-level data — which may be essential for the 
development of effective management strategies — are available, due to the lack of appropriate 
molecular markers. In order to contribute to the development of these strategies, we have developed 
the first set of microsatellite markers for P. moratoi. 

2. Results and Discussion 

2.1. Characterization of the Enriched Microsatellite Library 

A total of 384 clones were isolated and sequenced bidirectionally. The Codoncode Aligner 3.7.1 
software (CodonCode Corporation: Centerville, MA, USA) revealed a redundancy of 17% in the 
library. Of the unique clones selected for analysis in Microsatellite Repeats Finder [16], 176 (46%) had 
at least one microsatellite. A predominance of dinucleotide repeats (56%) was found in the motifs that 
make up the library. The CA W /GT W repeats (130 motifs identified) were the most numerous, followed 
by CTjv/GAat (48 motifs). This predominance of CA N /GT N repeats is typical of the eukaryote 
genome [17]. Considerable numbers of other types of motifs were also recorded, in particular the 
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dinucleotide AT W /TA W (36 motifs), the trinucleotides CATVGTAw (18 motifs), CTTVGAAyv (17), 
AATa/TTAjv and CTCVGAGyv (6 motifs each), and the tetranucleotides CTATVGATAw (13 motifs) 
and CATTyv/GTAAv (6 motifs). These data provide a baseline that will support the development of 
additional probes for the isolation of new microsatellites in P. moratoi. 

2.2. Development of Polymorphic Microsatellite Markers 

A total of 29 pairs of primers were designed and optimized successfully for the PCR amplification 
of microsatellite loci (Table 1). Of the loci analyzed, Pmoratoi\i\ presented several nonspecific 
amplifications even after optimization (with varying concentrations of magnesium chloride and 
different annealing temperatures) and was excluded. Six loci — Pmoratoi\i2, Pmoratoi\ii, Pmoratoi\iA, 
Pmoratoi\i9, Pmoratoi^lO, and Pmoratoi\i22 — were monomorphic. With the exception of 
Pmoratoi[i22, all these monomorphic loci represent interrupted or interrupted compound microsatellites 
characterized by a small number of repetitions, with predictably low polymorphism [18]. Twenty-two 
microsatellites were polymorphic (Table 2) in at least some populations {Pmoratoi\il , Pmoratoi\i%, 
Pmoratoi^ilO, Pmoratoi\i\\, Pmoratoi\i\A, Pmoratoi\i\l , Pmoratoi^ilS, and Pmoratoi\i2\). The 
Pmoratoi\i5 locus was not amplified in some populations, possibly due to local mutations in the primer 
annealing site. The identity analysis calculated using Cervus 3.0.3 [19] detected four pairs of 
specimens with identical genotypes (exact match), suggesting the recapture of the same animal during 
fieldwork. In these cases, the duplicate genotype was excluded from the analyses. 



Table 1. Microsatellites isolated in the present study with their respective primers and the 
optimal amplification conditions determined following visualization of the polymerase 
chain reaction (PCR) products in a polyacrylamide gel. 



Genbank 
Accession n' 


Locus 

> 


Repeat Motif 


Primer Sequence (5'^3') 


T A 


MgCl 2 (mIYl) 


JX441952 


Pmoratoi\il 


(TTTC) 9 


Forward: GGTGAACATCCTTTTCGTAGC 
Reverse: CACTCCTTCCCTAATCCAGTTT 


50 °C 


0.6 


JX441953 


Pmoratoi\i2 


(AC)4AT(AC) 7 
(AC) 4 


Forward: ACACATCGTTCTGCACTACACAC 
Reverse: GCTCCCTTGTCTTGCTGTCT 


63 °C 


1.0 


JX441954 


Pmoratoi\ii 


(TA) 8 CACA 
CAT(AC) 8 


Forward: CTAACCGTCCAATAGCCTGTGT 
Reverse: CCTCTTTCCCCTTGTGTGTCT 


63 °C 


0.6 


JX441955 


Pmoratoi\iA 


(AC) 6 G(CA) 5 


Forward: AAATGAGGTGGCTGTGCTAAAT 
Reverse: ATGCATTAGTGGTCATCACTGG 


60 °C 


3.5 


JX441956 


Pmoratoi\i5 


(CA) 8 


Forward: TATCTGTATTGCCTGCTCCACAC 
Reverse: CCTAGTGAGCTAAAAGTTGTGCTTGT 


68 °C 


3.5 


JX441957 


Pmoratoi\i6 


(ACAT) 4 
(AC) 15 


Forward: CTGCACCACCCCTTGAATAA 
Reverse: TGCACAGCAGGATCAATCTAAC 


46 °C 


0.8 


JX441958 


PmoratoiyO 


(AC) 8 


Forward: ACTTCCAGGTGCCATATCTTCA 
Reverse: AATTCTTGGTCTGCCATACTGTG 


51 °C 


1.0 


JX441959 


Pmoratoi\& 


(AC) 6 


Forward: GCGAATAATTGGAAAGCACAG 
Reverse: GCCTGAGCCAGAGTTGAATAGTA 


68 °C 


3.5 


JX441960 


Pmoratoi\i9 


(ATT) 4 ...(TAT) 4 


Forward: GATAATTGACCGTTTCCGTCAT 
Reverse: CATGGAACAAACTGAAGAGAACC 


63 °C 


4.0 
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Table 1. Cont. 



Genbank 
Accession n 


Locus 

> 


Repeat Motif 


Primer Sequence (5'— >3') 


T A 


MgCl 2 (mM) 


JX441961 


Pmoratoi\i\Q 


(TA) 4 (CA) 7 


Forward: CTAATAAAGTGGCCGGTGAGTG 
Reverse: ATAGGACTACATTGTGCCCTTGA 


50 °C 


0.8 


JX441962 


Pmoratoi\ill 


(CA) 8 


Forward: TCCAAAGTTCTAGCCTTGTTAG 
Reverse: CGCTACACATACCTTGAGAAA 


57 °C 


4.0 


JX441963 


Pmoratoi\i\2 


(ATCT) 7 
(CA) 5 ...(CA) 4 


Forward: CCTTCCCACCTTCCCTCTC 
Reverse: CGATCAACCTCCTCTTCTGTCTAC 


66 °C 


2.5 


JX441964 


Pmoratoi\il3 


(CA) 7 


Forward: CTGTTTGGACTGCGATTCTT 
Reverse: GCATTTGTGTGTGAGAGTGAA 


50 °C 


1.0 


JX441965 


Pmoratoi\i\A 


(ACAT) 8 


Forward: GTCAAATGAGGCGGCTGTG 
Reverse: GCCATTATTGCTTGTATTGCTTCAG 


63 °C 


1.5 


JX441966 


Pmoratoi\i\5 


(GATA) 12 
(CA)g 


Forward: CTTTAGGGCAGTCCAAGATTA 
Reverse: TGAAGGGGACACATTTTAAG 


50 °C 


1.5 


JX441967 


Pmoratoi\i\6 


(TCA) 8 


Forward: CTACACTAAAACGTCTCAATCAATG 
Reverse: ATGATGAAGAACTGGAGGAAGA 


66 °C 


2.5 


JX441968 


Pmoratoi\ill 


(CAC) V 


Forward: CCCAAAGAGTGCCAAGAAAATA 
Reverse: GGTAACAAACAACAAACCAGTATCAAC 


60 °C 


0.8 


JX441969 


Pmoratoi)y\% 


(CA) 7 


Forward: GTGTAATCCTGGGGTTCAGGTA 
Reverse: TCCCACCTTGGTCAGATATTGT 


57 °C 


1.0 


JX441970 


Pmoratoi\i\9 


(CA) 8 


Forward: TATAGTCCAGGCAGCCCCTTTA 
Reverse: GTCCGTGAGTGACGCAAAGT 


68 °C 


1.0 


JX441971 


Pmoratoi\i2Q 


(CA) 5 AG(CA) 6 


Forward: GATTCCCAGCAGAACATCAC 
Reverse: GGACTATGGAGCAATGAAAGAA 


63 °C 


0.8 


JX441972 


Pmoratoi\i2l 


(CA) 4 ...(CA) 4 ... 
(CA) 5 ...(CA) 7 


Forward: GGGGCACAGTGTATATGTCAGT 
Reverse: TTGAGCTGGTGAGGCAGTT 


66 °C 


3.0 


JX441973 


Pmoratoi\i22 


(TTTC) 17 


Forward: AAAATTCCGCTCAGTCATTA 
Reverse: ACTCCTTCCCTAATCCAGTTT 


42 °C 


4.0 


JX441974 


Pmoratoi\i23 


(TA) 4 (CA) U 


Forward: ACCTGGTCTAACCCTTTGGAAAT 
Reverse: CAGCGTTACCAGACATTTTATGTTC 


70 °C 


3.0 


JX441975 


Pmoratoi\i2A 


(AT) 7 


Forward: GCTATTTGTCTACCTATCTATCTTTCAT 
Reverse: CAATAAAACTCTGGACCTTGAAC 


40 °C 


4.0 


JX441976 


Pmoratoi\i25 


(ACT) n 


Forward: TCTAATGTCCACACTGCTACTACT 
Reverse: GCTAATGGCCGAGTTATTG 


70 °C 


4.0 


JX441977 


Pmoratoi\i26 


(CA) 7 


Forward: ATTTGGCTGTCTGACCTGTCTTA 
Reverse: CCCATATTAGTTCGGATCACAAG 


63 °C 


4.0 


JX441978 


Pmoratoi\i21 


(TCTA) 17 


Forward: CTCTATCTAACCCTTTCATA 
Reverse: AAGATGGATAGATGTGAGA 


57 °C 


2.0 


JX441979 


Pmoratoi)\2% 


(CA) 8 


Forward: GAAATGAGAGGCGTGAGAGAT 
Reverse: GCTGTCCGTCAATGGGTAT 


51 °C 


1.0 


JX441980 


Pmoratoi\i29 


(CA)i« 


Forward: GAGGAAAAGTCAAGGAACTAAATGTC 
Reverse: ACAGTCTTCTCAATCTGCATGTCT 


46 °C 


0.8 



T A : annealing temperature; MgCl 2 : magnesium chloride. 
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Table 2. Descriptive analysis of the genetic diversity in 22 polymorphic microsatellite loci 
obtained from five populations of P. moratoi. Significant deviations (p < 0.002) from 
Hardy- Weinberg Equilibrium (HWE) following the Bonferroni correction are indicated by 
an asterisk (*). Heterozyogosity and HWE were not estimated for the populations with 
small sample size (Avare and Lencois Paulista). 



Population 
Locus 




Sao Carlos 
(J! = 41) 




(n = 71\ 






I>I UlaS 






LP 

f n = 1> 
I" J ^ 


Total 




■''A 


H 0 


He 


/V 


ZJ 

''() 






ZJ 

n o 


it 

1 ' F 


J y a 


?V 
J' A 


S 


N A 


r 1 v 


Pfnorcitoi\iS 




0.29 


0.35 








-5 
J 


n ha 
u.uo 


ft ftf. 
u.uo 






201-239 


6 




Pfnovcitoi\i6) 


Q 

o 


0.70 


0.71 


C 

D 


u.oz 


U.Oo 


V 


u. jy 


n ^8 


z 


1 
1 


208-240 


10 


n 78 

U. /o 


Pwiorcitoi\il 


1 
1 


0.00 


0.00 


1 
1 


C\ C\C\ 

u.uu 


u.uu 


Z 


f\ /IA 

u.40 


ft A e 


i 
1 


1 
1 


185-187 


2 


a. 78 
U.zo 


Pfnorcitoi\L < & 


7 

Z 


0.17 


0.16 


1 
1 


u.uu 


u.uu 


z 


u.uz 


n 0.7 
u.uz 


1 
1 


1 
1 


209-211 


2 


ft ftf. 

u.uo 


Pmoratoi\i\Q 


n 
j 


0.30 


0.31 


1 
1 


u.uu 


U.UU 


-5 
J 


A IT 

U. 1Z 


n i 7 
u.lz 


1 
1 


1 
1 


109-113 


3 


ft 1 ^ 
U.l J 


Pmoratoi\i\ 1 


7 
z 


0.05 


0.05 


1 
1 


u.uu 


U.UU 




U. 1Z 




1 
1 


1 
1 


125-129 


3 


ft ftl 
U.U / 


Pmoratoi\i\2 


0 


0.80 


0.77 


Q 
0 


U.Dz 


U.o4 


1 1 
13 


U.oj 


U.6 / 


i 


Z 


144-196 


14 


ft 87 
U.6Z 


PmovatoiyL 1 3 


7 

z 


0.12 


0.16 


-5 
J 


u.^-z 


u. jy 


z 


U.U J 


ft fK 

U.U J 


z 


1 
1 


162-168 


4 


U. 


Pwiorcitoi\i\^ 


n 

1 


0.75 


0.76 


1 
1 


u.uu 


C\ C\C\ 
U.UU 


c 

J 


n ag 

U.Oo 


a. 7 1 


i 
1 


z 


190-218 


8 


n 7 1 


Pwiorcitoi\i 1 5 


i ft 
1 u 


0.67 


0.80 


o 

y 


\J.OJ 


n 7Q 
U. ly 


1 ? 

1Z 


u. / J 


ft 


z 


z 


177-241 


18 


n 84 


Pmoratoi\\,\6 


3 


0.62 


0.55 


3 


0.75 


A £ A 

0.64 


z 


A 1 A 

0.10 


0.10 


z 


1 


146-158 


4 


A /I O 

0.48 


Pmoratoi\ill 


3 


0.08 


0.08 


2 


0.04 


0.04 


1 


0.00 


0.00 


1 


1 


092-098 


3 


0.09 


Pmoratoi]y\% 


1 


0.00 


0.00 


1 


0.00 


0.00 


1 


0.00 


0.00 


2 


1 


163-165 


2 


0.02 


Pmoratoi]\,\9 


2 


0.50 


0.51 


2 


0.46 


0.49 


2 


0.36 


0.43 


2 


1 


167-169 


2 


0.40 


PmoratoiyClX 


1 


0.00 


0.00 


2 


0.40 


0.47 


2 


0.17 


0.16 


1 


1 


242-246 


3 


0.25 


Pmoratoi\iZi 


9 


0.72 


0.86 


4 


0.29 


0.61 


7 


0.76 


0.81 


3 


1 


225-251 


12 


0.87 


Pmoratoi\i2A 


3 


0.26 


0.66* 


2 


0.30 


0.47 


2 


0.10 


0.09 


1 


2 


148-174 


8 


0.72 


Pmoratoi\i25 


3 


0.60 


0.55 


3 


0.57 


0.58 


2 


0.54 


0.51 


1 


2 


246-252 


3 


0.53 


Pmoratoi\i26 


2 


0.13 


0.12 


2 


0.42 


0.38 


2 


0.41 


0.46 


1 


1 


119-133 


4 


0.38 


Pmoratoi)i21 


8 


0.45 


0.77* 


8 


0.87 


0.84 


9 


0.69 


0.85 


4 


2 


196-244 


12 


0.86 


Pmoratoi\i2S 


2 


0.45 


0.49 


3 


0.58 


0.58 


3 


0.61 


0.62 


1 


1 


195-203 


3 


0.48 


Pmoratoi\i29 


8 


0.90 


0.84 


7 


0.96 


0.77 


6 


0.47 i 


0.64 * 


3 


1 


154-198 


11 


0.80 



LP: Len9ois Paulista; 0.00: monomorphic locus; -: locus not amplified; N A : number of alleles; H 0 : observed 



heterozygosity; H E : expected heterozygosity; S: size range; PIC: polymorphic information content. 

The total number of alleles per locus (Na) varied between 2 and 18 (Table 2). Observed 
heterozygosity (Ho) ranged from 0.02 to 0.96, expected heterozygosity (H E ) from 0.02 to 0.87, and 
polymorphism information content (PIC) from 0.02 to 0.87. As might be expected from the 
relatively large number of markers developed for this study, evidence of linkage disequilibrium 
was found in three pairs of loci (Pmoratoi\i\0-Pmoratoi\i25 , Pmoratoi\i\2-Pmoratoi\i\5 , and 
Pmoratoi\i\5-Pmoratoi\i25) following Bonferroni correction (p < 0.002). Significant deviations from 
Hardy-Weinberg Equilibrium (HWE) were found in Pmoratoi[i24 and Pmoratoi\i27 from the Sao 
Carlos population and in Pmoratoi\x29 from Brotas, due to a deficit of heterozygotes. These deviations 
can be attributed to the presence of null alleles in these populations. The estimated null allele 
frequency for the Pmoratoi\i2A locus from Sao Carlos was 0.23, and that for Pmoratoi\i21 from this 
same population was 0.18. The estimated frequency for Pmoratoi\i29 from Brotas was 0.11. 
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In addition to the loci with deviations from HWE, null alleles were detected in the Pmoratoi\i\3 
(0.10), Pmoratoi\x23 (0.19), and Pmoratoi\x24 (0.10) loci from Bauru, and in PmoratoiyCll from Brotas 
(0.07). In all these cases, the evidence of the presence of null alleles was relatively weak and thus 
insufficient to confirm a significant departure from HWE following the Bonferroni correction. 
Micro-Checker 2.2.3 [20] did not detect small allele dominance, but found evidence of the presence 
of stutter bands in one locus, Pmoratoi\i\l> from Bauru. The identity analysis indicated that the 
combination of all the loci would permit the individual identification of each of the specimens. 

3. Experimental Section 

3.1. Construction of Enriched Micros atellite Genomic Library 

We constructed an enriched partial microsatellite genomic library using an approach based on the 
selective hybridization method of Kijas [21]. The library was constructed using DNA extracted from 
the muscle tissue of one specimen of P. moratoi using the procedure of Sambrook et al. [22] with 
modifications. Six micrograms of genomic DNA were digested with 50 units of Afa I (Invitrogen) and 
the fragments were then ligated to Rsa I linkers {Rsall: 5'-CTCTTGCTTACGCGTGGACTA-37 
Rsa25: 5 '-TAGTCC ACGCGTAAGC AAGAGC AC A-3 ') using 2 units of T4 DNA ligase (Promega). 
The fragments were then amplified by polymerase chain reaction (PCR) with a reduced number of 
cycles (20 cycles) using the primer Rsall. The PCR products were purified, denatured and hybridized 
with biotinylated microsatellite probes (GTg and CTg) at room temperature for 20 min. The hybrid 
mixtures containing microsatellites were then collected by streptavidin-coated magnetic beads 
(Promega). The selected fragments were amplified via PCR and the products were ligated into a 
pGEM-T easy cloning vector (Promega). Escherichia coli XL 1 -Blue cells (Stratagene) were 
transformed with recombinant plasmids by electroporation and grown overnight in solid Luria-Bertani 
agar medium containing ampicillin, IPTG and X-Gal. The positive colonies were selected and grown 
in liquid medium with 2YT HMFM containing ampicillin. After growing for 16 h, they were stored 
at -80 °C. 

3.2. Sequencing and Primer Design 

Of the total of 596 clones obtained, 384 were sequenced bidirectionally in an ABI Prism 3100 
automatic sequencer (Applied Biosystems: Foster City, CA, USA). The DNA sequences were exported 
into Codoncode Aligner 3.7.1 (CodonCode Corporation) which assembled the contigs and verified the 
redundancy of the library. The Bioedit program was used to check the quality of the sequences by 
chromatogram and to align them to form a consensus sequence. The repetitive elements were located 
using the Microsatellite Repeats Finder program [16]. After removal of the vector sequences, adapters, 
and restriction endonuclease sites by the Microsat software (version 1.0; CIRAD: Montpellier, France, 
2005), the primers were designed using Primer 3 [23]. 

3.3. Genotyping 

The polymorphic microsatellite markers were characterized by the amplification of the genomic 
DNA obtained from buccal epithelial cells (non-destructive method) following a modified version of 
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the procedure described by Pidancier et al. [24]. Samples were obtained from five remnant P. moratoi 
populations in the Brazilian state of Sao Paulo: 41 samples were collected in the municipality of 
Sao Carlos (22°01'00.5" S, 47°56'21.0" W), 41 in Brotas (22°12'53" S, 47°54'41" W), 27 in Bauru 
(22°20'48.46" S, 49°0'56" W), 3 in Avare (22°53.227' S, 48°56.803' W), and 1 in Lencois Paulista 
(22°49'13.17" S, 48°53'0.28" W). The PCRs were prepared in a final volume of 15 uL containing 10 ng 
of the DNA template, lx reaction buffer, 0.3 mM dNTP, 0.6-4.0 mM MgCl 2 (Table 1), 0.6 uM of each 
primer, and 1 unit of Taq polymerase (Invitrogen). The reactions were conducted following the same 
cycling conditions: 5 min at 94 °C followed by 41 cycles of 30 s at 94 °C, 1 min at the locus-specific 
annealing temperature (Table 1), and 1 min at 72 °C, followed by a final extension of 30 min at 72 °C 
to minimize stutter bands. The PCR products were analyzed in a Dual Dedicated Height Sequencing Kit 
(CBS Scientific) vertical electrophoresis system in 6% denaturing polyacrylamide gel and stained with 
silver nitrate [25]. Allele size was estimated by comparison with a 10 bp DNA ladder (Invitrogen) and 
using the GelAnalyzer 2010a software [26]. 

3.4. Characterization of Polymorphic Markers 

The levels of polymorphism of the microsatellites were evaluated as the number of alleles per 
locus (N A ), observed (Ho) and expected (H E ) heterozygosity calculated by Popgene 1.32 [27]. The 
polymorphism information content (PIC) was calculated with Cervus 3.0.3 [19], which was also used 
to conduct a test of individual discrimination (identity analysis). The Genepop 4.0.9 software [28] was 
used to detect an excess or deficiency of heterozygotes, linkage disequilibrium between pairs of loci, 
and deviations from the Hardy-Weinberg Equilibrium (HWE), for which significance levels were 
determined using the Markov chain algorithm [29], with 10,000 dememorization steps, 100 batches 
and 5000 iterations per batch. All significance levels were adjusted by the sequential Bonferroni 
correction for multiple tests [30]. Micro-Checker 2.2.3 software [20] was used to identify genotyping 
errors and verify the presence of null alleles using the Brookfield method [31]. 

4. Conclusions 

These are the first microsatellite markers developed for Morato's digger toad, and in fact, the first 
for any member of the genus Proceratophrys. These markers constitute a powerful tool for the study of 
P. moratoi populations, allowing the identification of untagged individuals and providing a database 
for the development of kinship studies for future ex situ conservation programs. It will also be possible 
to analyze inbreeding, genetic diversity and structure, and gene flow in natural populations, which will 
be vital for the development of effective in situ conservation measures. 
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